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Pairs potential-density in terms of elementary functions that represents flat rings structures are 
presented. We study structures representing one or several concentric flat rings. Also disks sur- 
rounded by concentric flat rings are exhibited. The stability of concentrically circular orbits of 
particles moving on a flat ring structure is analyzed for radial perturbations. 



o 
o 

(N 
C 



6 ' 



> 

O 

oo 

o 
i> 
o 



X 



I. INTRODUCTION 

Flat rings systems are a common feature of the four 
giant planets of the Solar System. Also ringy structures 
are present in several observed galaxies and nebulae, e.g., 
in the Hourglass nebula we have as an outstanding char- 
acteristic: the presence of a very net and large isolated 
ring. The exact gravitational potential of a ring of zero 
thickness and constant linear density is given in its exact 
form by an elliptic integral that is seldom used for practi- 
cal purposes. This potential is usually approximated by 
truncated series of spherical harmonics, i.e., a multipolar 
expansion. As far as we know there not simple expres- 
sions for the exact gravitational potential of a flat ring of 
any surface density. 

The potential of a ring enclosing a disk can be ob- 
tained by a process of c omp lcxification of the potential 
of a punctual mass 

Q> & 11 • This potential con be 
used to built a family of similar structures @. On the 
other hand simple pairs potential-density for thin disks 
are known. The simpler is the Plummer-Kuzmin disk 
[ fPH and [||] that represents a simple model of galactic 
disk with a concentration of mass in its center and den- 
sity that decays as 1/r 3 on the plane of the disk, see 
also Q . This structure has no boundary even though for 
practical purposes one can put a cutoff radius wherein 
the main part of the mass is inside, say 98% of the mass. 
Another simpler models of disks are the [l(| disks. These 
disks have a mass concentration on their centers and fi- 
nite radius. The [l2| inversion theorem can be used to in- 
vert the Morgan and Morgan disks to produce an infinite 
disk with a central hole of the same radius of the original 
disk. We can also put a cutoff in these inverted disks. 
Therefore the inverted Morgan and Morgan disks can be 
considered as representing a flat rings. In the context 
Einstein theory of gravitation one of this inverted rings 
was used to study the gravity of a disk with a central 
black hole, 0]- 

The purpose of this article is to use the family of Mor- 
gan and Morgan disks to construct flat ring like struc- 



tures. We shall use two different approaches, the first 
consists in the superposition of finite disks of different 
densities. By using this method we can construct struc- 
tures that represent one or several concentric flat rings 
and families of disks surrounded by flat rings. These 
structures have a finite outer radius. The second method 
is the Lord Kelvin inversion method that will be used to 
construct the corresponding inverse structures obtained 
using the first method. In this case the structures extend 
to infinite, but as before one can also put a cutoff. 

As a very first test of stability for these structures we 
study the linear stability of particles moving in circular 
orbits inside them. In other words, we assume the very 
naive model that the structures are build from particles 
moving only in concentric circles. The stability of circular 
orbits in several axially symmetric systems is studied in 
both Newton and Einstein theories of gravitation in Q . 

This article is divided as follows, in Sec. II we give 
a quick revision of the pair potential-density associated 
to the Morgan and Morgan disks. In Sec. Ill we study 
several classes of superpositions that give rise to pairs 
potential-density of structures that can be considered 
as representing one or several concentric flat rings. In 
the next section, Sec. IV, we construct pairs potential- 
density associated to structures representing a central 
disk surrounded by one or several concentric flat rings. 
In the following section, Sec. V, we study the Kelvin in- 
version of the previous studied structures. In the penulti- 
mate section, Sec VI, we study the stability under linear 
radial perturbations of four of the ring structures already 
presented. Finally, in Sec VII, we summarize and discuss 
some of the previous results. 



II. MORGAN AND MORGAN DISKS 

The [l(| disks are obtained by solving the Laplace 
equation in the natural coordinates to represent the grav- 
itational potential of a disk like structure, i.e., prolate 
coordinates (v, (i, if) that are related to the usual cylin- 
drical coordinates (r, z,tp) by 
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av = — ImT?., a[L = RelZ, 1Z = \Jr 2 + (z — ia) 2 , (1) 



2 



where a is a positive constant. From the previous equa- 
tions we find 

r = a^Jv 2 + 1)(1 - fi 2 ), z = av/j,. (2) 

Note that on the plane z — we have v = and fi = 
y/l — p 2 , with p = r/a. 

From the Laplace equation in prolate coordinates one 
can find, 0], the potential of a disk parallel to the plane 
z = and centered on the origin of the coordinate system, 

V 2m = (1 - 2m)GM 2m (~l) m P 2m (0)P 2m ( f i)q 2m (v)/a, 

where M 2m is the mass of the disk with potential V 2m and 
a its radius. P n {x) are the usual Legendre polynomials 
and q n (x) are related to the usual Legendre functions, 
Q n (x), by q n (x) — i n+1 Q n (ix). We recall the identity, 

(-l) n P 2n (0) = 1 -3-5 ... (2n- l)/(2-4-6 ... 2n), 

= (2n-l)!!/(2n)H, (4) 

that will frequently used. 

Associated to the potential (J3J) we have the surface 
density, @, 

o 2m = (2m + l)Af 2m F 2m (l -p 2 )/(2vraVl -P 2 )- (5) 

[HI considered the following superposition of the above 
mentioned solutions of Laplace equation, 

rn—n 

V (n) - 51 A 2 m ,2n^ 2m , (6) 
m=Q 

with 

(4m+l)2 2m (2n)! 
2m ' 2n (n + m)!(n-m)!(2n + 2m + l)!' 1 ' 

This superposition represents the potential of disks with 
the simple associated surface density, 

a^(r) = (2n+ l)Af (,l) (l -p 2 )™^/^™ 2 ), (8) 

where we have taken M 2n = M^> for all n. 

The potential on the plane z — can be obtained using 
the algorithm, [13], 

= (2n+l)(-l) n P 2 „(0) / [1-p 2 cos 2 6] n d9 

a Jo 

(9) 

The n = disk is singular on its rim and using (TTJ) can be 
interpreted as the potential associated to a complexified 
bar of constant density, @. The subsequent members 
of the family represents disks of finite density with their 
maximum on their centers and zero density on their rims. 

III. FLAT RINGS 

Let us consider disks of the same radius a and decreas- 
ing mass, 

M (n) = 2Tra c a 2 /(2n+ 1), (10) 
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FIG. 1: The surface density of the first three members of a 
family disks with a central hole, {o~[™ }, i.e., flat rings, the 
size of the central hole is bigger for larger n. 

where a c is a constant and will be taken equal for all 
disks of the Morgan and Morgan family, n = 1,2,3,.... 
For these family of disks we have that the corresponding 
surface density is 

<J {n) =a c (l-p 2 )'(l-p 2 )™- 1 . (11) 
Now let us consider the following superposition, 

n 

°lr = ]rC£(-l) n - fc a(" +1 - fc >, (12) 

fe=0 

= a c (l-p 2 )h 2 \ (13) 

where C£ — n\/[(n— k)lk\]. We have that all these super- 
positions give disks of radius a with zero density on their 
centers, i.e., disks with a hole in their centers, in other 
words flat rings. The densities associated to the three 
first members of this family of flat rings are presented in 
Fig. [TJ We see a disk with a hole in the center with a 
residual density that is small for larger n. The mass of 
the three first flat rings are 7rcr c a 2 /16, 7rcr c a 2 /32, and 
57TCT c a 2 /256, respectively. The potentials associated to 
these flat rings can be found using a superposition with 
the same coefficients as the ones used for the densities, 
e.g., the potential associated to <j\ r is V\ r = — V^ 2 \ 
etc. 

Now let us consider, instead a superposition of disks, 
a superposition of flat rings, 

(n) (n) ,2 («) , iA ( n ) /i a\ 

where b is a constant such that b > a. This superposition 
put a gap at p — 1/6 (r = a/b) on the flat one ring given 
us a family of two concentric flat rings. In Fig. [2] we 
present and cr^J , with b 2 — 2, as before we see that 
the larger the n the larger is the size of the central hole. 
In Fig. OH we show a$ /a c for b 2 = 2, 1.7, 1.4. 

We can construct several families of two rings using 
one rings in such a way that the final density be 

4 n r ^=a c (l-p 2 )Hl-bVf m P 2n (15) 
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FIG. 2: The surface density of the first two members of a FIG. 4: The surface density of the first member of the family 



family disks with a center hole and a gap, a^} and a^.' , i.e., of disk-one-rings, cr^ , for b z = 3, 2.7, 2, 4, 2.1. For b small we 



(i) 



flat double rings, the size of the central hole is bigger for larger 
n. We take b 2 = 2. 



have a larger central disk. 
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FIG. 3: The surface density of the first member of the family 
of two-rings, <r£?, for 6 2 = 2, 1.7, 1.4. 



Of course this kind of superposition can be used to have 
families of solutions representing n concentric flat rings 
with density, 

at^ = a c (l-p 2 )i(l - bf 1)P 2 ) 2mi a - bf 2) p 2 ) m2 x 
.■■x(l-bl_ 1)P 2 ) 2m ^p 21 , (16) 

with m = (mi, m„_i) S N"^ 1 . The positions of the 
gaps are given by pi = 1/6,:, with i = l...n— 1, and 6,; > a. 

IV. DISKS WITH FLAT RINGS 

By using the same type of superpositions presented in 
the previous section we can construct the pair potential- 
density for disks surrounded by flat rings. Let us consider 
models of a disk surrounded by one concentric flat ring. 
The density 

= ^-2^%+^, (17) 
= a c (l-p 2 )Hl-b 2 p 2 ) 2 , (18) 



represents a disk of radius r = a/b with a flat ring be- 
tween r — alb and r — a. In Fig. 0]we present cr^ r /<j c 
for b 2 = 3, 2.7, 2, 4, 2.1 we see that the position of the gap, 
as well as, the maximum of the fiat ring density depends 
on b. 

As before, we can build a family of disks with sur- 
rounded by one ring with surface density, 

^ r ^a c (l-p 2 )Hl-b 2 p 2 r, (19) 

and the general case of a disk surrounded n concentric 
flat rings belonging to the above mentioned family, 

••.x(l-^/) 2 ™"{20) 
where now rh = (mi, m n ) S N n . 

V. KELVIN INVERTED FLAT RINGS 

To construct new family of flat rings we can use the 
inversion theorem of [12j that in cylindrical coordinates 
tells us that: if the pair potential-density V(r, z), cr(r, z) 
is a solution of Poisson equation the pair, 

V(r, z) = V [a 2 r/(r 2 + z 2 ),a 2 z/{r 2 + z 2 )} , (21) 
a{r, z) = {a/rfa [a 2 r/{r 2 + z 2 ), a 2 z/(r 2 + z 2 )](22) 

is also a solution of the same equation, see for instance, 
Q . The inversion of the Morgan and Morgan family of 
disks with the density §E§ gives us, 

<7 (n) (r) = {2n+l)M^/(2na 2 q 3 )(l-q 2 ) n -i, (23) 

where q = a/r and n > 1. Note that that the Kelvin 
inverted Morgan and Morgan disks can be considered 
as representing flat rings, we have a disk with a hole 
of radius a. The density of these fiat rings decays in a 
similar way as the Plummer-Kuzmin disks, i.e, as 1/r , 
not too fast. 
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FIG. 5: The surface density of the first three member of the 
family of inverse one-rings, {rrj"'}. We have a central hole of 
radius r — a. Even though the rings extent to infinity we see 
a clear cutoff. 



The mass of the Kelvin inverted Morgan and Morgan 
disks is given by, 

M (n) = (2n + l)!!(2n)!!M (Tl) /(2n)!!. (24) 

The Kelvin inverted one-rings (|13[) gives other family of 
one-rings with density 

a^=a c (l- q ^q 2n+3 - (25) 

Note that the density decays as <r 1 ™' ) ~ i/ r 2 ™+ 3 for this 
family of flat rings. In Fig. [5] we plot the inverted den- 
sity of the first three members of the one-flat-ring family. 
Note the central hole of radius r — a. Due to the fast 
decay rate with the distance we have that on can put a 
clear cutoff and to consider this flat rings as finite, Eq. 
(|2l))) . We can also invert the families of two and sev- 
eral flat rings discussed in Sec. 3, we will also have new 
families of rings that decay with the distance very fast. 



VI. STABILITY OF CIRCULAR ORBITS 

In this section we study the stability under radial per- 
turbations of circular orbits concentric to the flat rings. 
By assuming that flat rings are made of particles moving 
along circular orbits the study of the stability of these 
orbits can be considered as an order zero test of stability 
of the flat rings. In this case the collective behavior of 
the particles of the ring is not taken into account. The 
epicyclic frequency of the perturbed orbit is, [H, 

\ dp 2 p dp J 

The epicyclic frequencies of rings with densities 
a vr ' ano - their associated Kelvin inverted flat single 



rings with densities a^J , a^} are: 

«£> = 7T [a 3 a c (27p 2 - 8)/8] 1 , (0 < p < I), (27) 

K ir = \ [a 3 <Tc(105p 4 -8Ip 2 + 20)]^ , (0<p<l$8) 

/#? = 3^[a 3 a c /8)]^, (p > 1), (29) 

k£> = ^[3a 3 a c (-I8p 2 + 35)/32p, (p > I). (30) 

The one-ring with density a[]} is not completely sta- 
ble, we have stability only for r > (2v2/3)a. The 
flat rings with densities, af^ \ and a[]} are stable in 
all their extensions. We have that only the inner part, 
I < (r/a) < ^/35/I8, of the flat ring with density a^) is 
stable. 



VII. DISCUSSION 

We presented several families of pairs potential-density 
that can be used to represent one or several concentric 
rings, also a disk surrounded by flat rings. We have 
shown in some detail the superposition of surface den- 
sities to build the families indicating that the same su- 
perposition is valid to build the associated potential. 

We have considered only the density as representing 
the physics inside the different structures studied along 
this work, in principle, one can find the other thermody- 
namics variables by solving the Fokker-Plank equation. 
In a recent work we presented a method for solving this 
equation for a disk, [l3|, that with simple modifications 
can be used to find a distribution function for these struc- 
tures. 

With our very first approach to the stability of fiat 
rings presented in the penultimate section, we find sta- 
ble configurations, as well as, partially stable ones. The 
structure of one-ring standing alone is not very realist. 
For instance in the case of planetary flat rings we have a 
planet located in center of the ring that will contribute 
to stabilize the structure. To be more specific let us con- 
sider a planet with mass M surrounded by a flat ring. 
Let k 2 . be the square of the epicyclic frequency associ- 
ated to the ring alone. Now the square of the epicyclic 
frequency of a particle orbiting the system planet-ring 
is = ^ r 

+ GMa 3 /p 3 . Since the contribution of the 
planet to the epicyclic frequency of the particles of the 
disk is a positive quantity we have that the planet tends 
to stabilize the rings. To perform a study of the stability 
of the structures above described taken into account the 
collective behavior of the particles we need to know the 
fluid dynamics variables that in principle can be found 
solving the Fokker-Planck equation using the method al- 
ready mentioned. 
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